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THE INFLUENCE OF AN AMBIENT MAGNETIC FIELD ON RELATIVISTIC COLLISIONLESS PLASMA SHOCKS 
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ABSTRACT 

Plasma outflows from gamma-ray bursts, supernovae, and relativistic jets, in general, interact with the sur- 
rounding medium through collisionless shocks. The microphysics of such shocks are still poorly understood, 
which, potentially, can introduce uncertainties in the interpretation of observations. It is now well established 
that the Weibel two- stream instability is capable of generating strong electromagnetic fields in the transition re- 
gion between the jet and the ambient plasma. However, the parameter space of collisionless shocks is vast and 
still remains unexplored. In this Letter, we focus on how an ambient magnetic field affects the evolution of the 
electron Weibel instability and the associated shock. Using a particle-in-cell code, we have performed three- 
dimensional numerical experiments on such shocks. We compare simulations in which a jet is injected into 
an unmagnetized plasma with simulations in which the jet is injected into a plasma with an ambient magnetic 
field both parallel and perpendicular to the jet flow. We find that there exists a threshold of the magnetic field 
strength below which the Weibel two- stream instability dominates, and we note that the interstellar medium 
magnetic field strength lies well below this value. In the case of a strong magnetic field parallel to the jet, the 
Weibel instability is quenched. In the strong perpendicular case, ambient and jet electrons are strongly accel- 
erated because of the charge separation between deflected jet electrons and less deflected jet ions. Also, the 
electromagnetic topologies become highly non-linear and complex with the appearance of anti-parallel field 
configurations. 

Subject headings: acceleration of particles, gamma rays: bursts, shock waves, instabilities, magnetic fields, 
plasmas 



1. INTRODUCTION 

The collisionless plasma condition applies to many astro- 
physical scenarios, including, the outflow from gamma-ray 
bursts (GRBs), active galactic nuclei, and relativistic jets in 
general. The complexity of kinetic effects and instabilities 
makes it difficult to understand the nature of collisionless 
shocks. Only recently, the increase in available computa- 
tional power has made it possible to investigate the full three- 
dimensional dynamics of collisionless shocks. 

In the context of GRB afterglows, observations indicate 
that shock-compressed magnetic field from the interstellar 
medium (ISM) is several orders of magnitude too weak to 
match observations. Particle-in-cell (PIC) simulations have 
revealed that the Weibel two- stream instability is capable of 
generating the required electromagnetic fi eld strength of the 
order o f percents of e quipartiti on value jFre deriksen et al. 
(2004); Kazi mura et ajl (119981): iMedvedev & Loebl (11999): 
Nishikaw aet alJ (120031 l200a »: ISilva et atJ (I2003M Further- 
more, PIC simulations have shown that in situ non ther- 
mal particle acceleration takes place in the shock transi- 
tion region (Ei ededal et al. 2004; Hoshino & Shimada 2002; 
ISaito & S akai 2003|K Three-dime nsional simulation s using 
~ 10 7 electron-positron pairs by Sakai and Matsuo (2004) 
showed how complex magnetic topologies are formed when 
injecting a mildly relativistic jet into a force-free magnetic 
field with both parallel and perpendicular compon ents. With 
a two-dimensional analysis, Saito & Sakai (2003) found that 
an ambient parallel magnetic field can quench the two- stream 
instability in the weakly relativistic case. In this Letter, we 
use three-dimensional PIC experiments to investigate how 
the two- stream instability is affected by the presence of an 
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ambient magnetic field. Using up to ~ 10 9 particles and 
125 x 125 x 1200 grid zones, we investigate the development 
of complex magnetic topologies when injecting a fully rela- 
tivistic jet (bulk Lorentz factor F = 5) into an ambient magne- 
tized plasma. Using varying field strengths, we focus on the 
case of a transverse magnetic field and compare it with the 
case of a parallel magnetic field. The simulations are mainly 
concerned with the electron dynamics since processes involv- 
ing the heavier ions evolve on much longer timescales. 

2. THE NUMERICAL EXPERIMENTS 



We use the PIC code described by Frederikse n et alJ (120041). 
The code works from first principles and evolves the equa- 
tion of motion for the particles together with Maxwell's equa- 
tions. In the simulation experiments, we inject an electron- 
proton plasma (a jet) into an ambient plasma (the ISM) ini- 
tially at rest (Fig. [Q. The jet is moving with a relativis- 
tic velocity of 0.98c along the z-direction corresponding to 
Lorentz factor jj et = 5. The ion-to-electron mass ratio is set 
to mi/m Q = 20. The jet plasma and the ambient plasma have 
the same density, n, and the corresponding electron plasma 
rest-frame frequency uj pe = [ne 2 /(m e eo)] = 0.035 A" 1 (e is the 
unit charge, eo the vacuum permittivity, and A t the simula- 
tion unit time). We choose this low value in order to properly 
resolve the microphysics. Initially, the interface between the 
ambient and the injected plasma is located at z = 3A e , where 
X e is the electron skin depth defined as X e = c/uj pe = 28.6A X 
(c is the speed of light, and A x the grid size). The time 
step and grid size obey the Courant-Friedrichs-Levy condi- 
tion A t = 0.5A x /c. Both plasma populations are, in their re- 
spective rest frames, Maxwellian distributed with a thermal 
electron velocity v t h — 0.03c. This temperature allow us to 
numerically resolve the plasma Debye length with approxi- 
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Fig. 1. — Schematic example of the simulation setup. Jet plasma is ho- 
mogeneously and continuously injected in the z-direction throughout the x—y 
plane at z = 0. Inside, the box is populated by a plasma population, initially 
at rest. In this specific example, an ambient magnetic field is set up in the 
x-direction (perpendicular case). 



We consider three different ambient magnetic configura- 
tions: no magnetic field, a magnetic field parallel to the flow, 
and a magnetic field perpendicular to the flow. The magnetic 
field is initially setup to be homogeneous and at rest in the am- 
bient plasma. The experiments are carried out with 1 billion 
particles inside (125 x 125 x 1200) grid zones. In terms of 
electron skin depths, this corresponds to (4.4,4.4,42)A e . The 
boundary conditions are periodic in the direction transverse 
to the jet flow (x, y). In the parallel direction, jet particles are 
continuously injected at the leftmost boundary (z = 0). At the 
leftmost and rightmost z boundary, electromagnetic waves are 
absorbed, and we allow particles to escape in order to avoid 
unphysical feedback. The total energy throughout the simula- 
tions is conserved with an error less than 1%. 

3. RESULTS 

Initially, we run simulations with no ambient magnetic 
field and observe the growth of the Weibel two- stream 
instabi li ty also found i n prev io us work (Frederiksen et al. 
(2004); Kazi mura et all (119981) : iMedvedev • & I^ejT Tl999): 
Nishika waetal.1 (120031 12005b : ISilvaetatJ (|2003l) V The 
Weibel two- stream instability works when magnetic pertur- 
bations transverse to the flow collect streaming particles into 
current bundles that in turn amplify the magnetic perturba- 
tions. In the non linear stage, we observe how current fila- 
ments merge into increasingly larger patterns. The electro- 
magnetic energy grows to e# ~ 1%, where e# describes the 
amount of total injected kinetic energy that is converted to 
magnetic energy. 

3.1. Parallel Magnetic Field 

In the presence of a strong magnetic field component par- 
allel to the flow, particles are not able to collect into bundles 
since transverse velocity components are deflected. We have 
performed five runs with parallel magnetic fields correspond- 
ing, respectively, to uj pe /uj ce =40, 20, 10, 5, and 1 while keep- 
ing uo pe constant; uo ce = eB/(jj et m e ) is the jet electron gyrofre- 
quency. The resulting field generation efficiency can be seen 
in Fig. |2|at t = 2\u~ l e where the jet front has reached z = 23X e . 
In the case of uo pe /uo ce = 40, the Weibel instability overcomes 
the parallel field, and although initially slightly suppressed, 
it eventually evolves as in the case of no ambient magnetic 
field. Increasing the magnetic field to uo pe /uo ce = 1 effectively 
suppresses the instability. Thus, for an ISM strength mag- 
netic field (ujpe/ujce ~ 1500) parallel to the plasma flow, the 



Weibel two- stream instability will evolve unhindered, and the 
generated field will exceed the ISM field. We find from the 
simulations that it would take a milligauss strength parallel 
magnetic field to effectively quench the instability for a 7 = 5 
jet expanding in an environment with density like the ISM. 




z (electron skindepth) 

Fig. 2. — Growth of the Weibel two-stream instability for different 
strengths of the parallel ambient magnetic field at time t = 2\uj~\. Here we 
measure the effectiveness of the field generation as the average transverse 
magnetic field amplitude as a function of z. The solid line corresponds to 
w pe/ w C e = 40, the dotted line to 20, the dashed line to 10, the dot-dashed line 
to 5, and the triple-dot-dashed line to 1. The case of no ambient magnetic 
field is very similar to that of uo P e/oo C e = 40. The magnetic amplitude is in 
arbitrary units 

The left panel of Fig. [3] show the resulting electron mo- 
mentum distribution function for different values of uj pe /ujce 
at t = 30cj~]. Since the presence of a strong parallel mag- 
netic field suppresses the generation of a transverse magnetic 
field, there exists no mechanism that can heat the electrons 
and transfer momentum between the two electron popula- 
tions. Thus the jet plasma propagates unperturbed. Where 
there is no parallel magnetic field or only a weak magnetic 
field (pupe/uoce = 1500), we observe how the jet and ambient 
plasma is heated and how momentum is transferred between 
the two populations. 

3.2. Perpendicular Magnetic Field 

We have performed experiments with an ambient mag- 
netic field perpendicular to the jet flow (Fig. [IJ with 
field strengths corresponding to 6^^/^=1500, 40, 20 and 
5. By including the displacement current, one can derive 
the relativistic Alfven speed = c~ 2 + (y^n-rehyi ~ c /[\ + 

(CVM,)W^)77^] 1/2 , where v— 1 = B/[fi n(mi + 
)] 1 / 2 is the non-relativistic counterpart. From this we cal- 
culate the corresponding relativistic Alfven Mach numbers 
JjetVjet/vA = 6572, 175, 88, and 22. 

Again, the uj pe /uj ce = 1500 run has been chosen because it 
resembles the typical density and microgauss magnetic field 
strength of the ISM. We find that the magnetic field generated 
by the two- stream instability dominates the ambient magnetic 
field, and the result resembles that of no ambient magnetic 
field. Furthermore, as expected in both the parallel and per- 
pendicular cases, the electron momentum distributions (Fig. 
are very similar, except for a weak merging between the 
ambient and jet electrons in the perpendicular case. 

In the run with uj pe /uj ce = 20, the result differs substantially 
from the previous cases. With reference to Fig. |4j we de- 
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Fig. 3. — Normalized electron momentum distribution functions at time 30uj p ] . The left panel is for runs with an ambient magnetic field parallel to the injected 
plasma with uj pe /uj ce = 1 (solid line), 20 (dotted line), 40 (dashed line), and 1500 (dot-dashed line). The right panel is for runs in which the initial magnetic field 
is perpendicular to the inflow and uj pe /uj ce = 5 (solid line), 20 (dotted line), 40 (dashed line) and 1500 (dot-dashed line). The vertical line shows the injected 
momentum 7 = 5. The distribution functions are for electrons with z > 15A e . 
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Fig. 4. — Snapshot at t = 16lj~1 of the highly complicated magnetic field 
topology when a jet is injected into a plasma with an ambient magnetic field 
transverse to the jet flow. The bottom panel shows magnetic field lines in a 
subsection of the computational box (from z = 9X e to z = I5\ e . The top panel 
shows a schematic explanation in the x—z plane: Jet electrons are bent by the 
ambient magnetic field (region A). Due to the Weibel instability, the electrons 
bundle into current beams (region C) that in turn reverse the field topology 
(region B). This will eventually bend the jet beam in the opposite direction. 
[See the electronic edition of the Journal for a color version of this figure.] 



scribe the different stages of the evolution: Initially, the in- 
jected particles are deflected by the ambient magnetic field. 



The magnetic field is piled up behind the jet front, and the en- 
hanced magnetic fields bend jet electron trajectories sharply. 
This has two implications. First, the ions, being more mas- 
sive, will penetrate deeper than the deflected electrons. This 
creates a charge separation near the jet head that effectively 
accelerates both ambient and injected electrons behind the ion 
jet front as shown in Fig. |5] Second, the deflected electrons 
eventually become subject to the Weibel two- stream instabil- 
ity. This forms electron current channels at some angle to the 
initial direction of injection, as shown in the upper panel of 
Fig. |4| Around these current channels, magnetic loops are 
induced (Fig. |4l region C). Magnetic islands are formed and 
the ambient magnetic field is reversed behind the loops (re- 
gion B). In this region, we find acceleration of electrons in 
the x-direction. The activity in this region has similarities to 
reconnection, but it is beyond the scope of this Letter to in- 
vestigate this topic. 
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Fig. 5. — Electrons v z j plotted against their position in the shock at time 
equal 30uj~1 for the run with an ambient magnetic field transverse to the 
jet flow corresponding to uj pe /u)ce = 20. The jet electrons (black dots) are 
injected at z = with 7 = 5. The ambient electrons (gray dots) are initially at 
rest but are strongly accelerated by the jet. 
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In other regions, the ambient magnetic field is strongly 
compressed, and this amplifies the field strength up to 5 times 
the initial value. As a result, parts of the jet electrons are actu- 
ally reversed in their direction. This can be seen from Fig. 
which shows a phase- space plot of both ambient and jet elec- 
trons at t = 30uj~1 . We see several interesting features here. In 
the region z = (15 - 21) A e , we observe how ambient electrons 
are swept up by the jet. Behind the jet front, both ambient 
and jet electrons are strongly accelerated since the jet ions, 
being heavier, take a straighter path than the jet electrons, and 
this creates a strong charge separation. The excess of positive 
charge in the very front of the jet head is very persistent and 
hard to shield since the jet ions are moving close to the speed 
of light. Thus, there is a continuous transfer of z-momentum 
from the jet ions to the electrons. In the case of the perpen- 
dicular ambient field, more violent processes take place than 
in the case of the parallel field, which can be seen in Fig. |3] 
(right panel). Here we see that mixing of the two plasma pop- 
ulations is much more effective for the perpendicular case. 
However, the spectrum of the electrons' momentum is highly 
nonthermal, with strong acceleration of both jet and ambient 
electrons. The cutoff in electron acceleration depends on the 
magnetic field strength. The maxima (7V n « 10) at z = 20 X e 
in Fig corresponds to the cutoff shown by the dotted line in 
the right panel in Fig. |3] It should be noted that the current 
channels that are caused by the bent jet electron trajectories 
in the early time, as shown in Fig. 0] are also seen in Fig. 
The first current channels have moved to around z = 20 X e . At 
z = 15 X e , a second current channel is created by the deflected 
jet electrons. This periodic phenomenon involves the ions in a 
highly nonlinear process but is beyond the scope of this Letter 
and will be explained in a subsequent paper. 

4. CONCLUSIONS 

Using a three-dimensional relativistic particle-in-cell code, 
we have investigated how an ambient magnetic field af- 
fects the dynamics of a relativistic jet in the colli- 
sionless shock region. We have examined how the 
different ambient magnetic topology and strength af- 
fect the growth of the electron Weibel two- stream in- 
stability and the associated electron acceleration. This 
instability is an important mechanism in collisionless 
shocks. It facilitate s momentum transfer between colliding 
plasma populations ( Frederikse n et al.ll2004 iK azimura et al l 
1998; Medvedev & Loeblh 999UNishikawa et al .1120031 12005 : 
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ISilva et at1l2003l) and can accelerate electrons to nonthermal 
distributions ( Hede dal et al.l2004t iHoshino & Shimada[ 2002 ; 
Saito & Sakai 2003). Collisionless shocks are found in the 
interface between relativistic outflows (e.g., from gamma-ray 
bursts and active galactic nuclei) and the surrounding medium 
(e.g., the ISM). 

We find substantial differences between the cases of am- 
bient magnetic fields transverse and parallel to the jet flow. 
However, common for both cases is that it takes an ambient 
magnetic field strength much stronger than the strength of the 
magnetic field typically found in the ISM to effectively sup- 
press the Weibel two-stream instability. In the case of a par- 
allel magnetic field, uo pe /uo ce must be smaller than 5 to effec- 
tively suppress the instability. This res ult is in good ag r eemen t 
with two-dimensional simulations by Saito & Sakai (2003), 
and thus this limit seems independent of jj et . For a typical 
ISM density of 10 6 m~ 3 , this corresponds to a milligauss mag- 
netic field. We emphasize the importance of uo pe /uo ce as an im- 
portan t parameter for collisionless shocks, as was also pointed 
out by Shimada & Hoshino ( 200^). 

In the case of perpendicular injection, the dynamics are dif- 
ferent from the parallel injection. Here, the electrons are de- 
flected by the magnetic field, and this creates a charge sep- 
aration from the less deflected ions. The charge separation 
drags the ambient and jet electrons, and consequently they 
are strongly accelerated along the z-direction. Furthermore, 
due to the Weibel instability, current channels are generated 
around the ambient magnetic field, which is curled and lo- 
cally reversed. 

These simulations provide insights into the complex dy- 
namics of relativistic jets. Further investigations are required 
to understand the detailed physics involved. Larger simula- 
tions (above 10 9 particles) with longer boxes are needed to 
cover the instability domain of the ions, to investigate the full 
evolution of the complicated dynamics, and to resolve the 
whole shock ramp. 
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